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Discrete Dubrovin Equations and 
Separation of Variables for Discrete Systems 

F.W. Nijhoff 

Department of Applied Mathematical Studies 
The University of Leeds, Leeds LS2 9JT, UK 

Abstract 

A universal system of difference equations associated with a hyperelliptic curve is 
derived constituting the discrete analogue of the Dubrovin equations arising in the theory 
of finite-gap integration. The parametrisation of the solutions in terms of Abelian functions 
of Kleinian type (i.e. the higher-genus analogues of the Weierstrass elliptic functions) is 
discussed as well as the connections with the method of separation of variables. 



1 Introduction 

The method of finite-gap integration has been developed in the late seventies to deal with 
the periodic solutions of integrable nonlinear evolution equations, cf. e.g. ffl, 0, M, and the 
recent monograph Q. The Dubrovin equations arise in the theory of finite-gap integration, 
namely as the equations governing the dynamics of the so-called auxiliary spectrum or 
equivalently of the poles of the Baker- Akhiezer function. In the case of integrable models 
associated with a hyperelliptic curve they are typically given by a first-order system of 
coupled ordinary differential equations the form, M, 



in which R([i) denotes the discriminant of the hyperelliptic curve in question. Variations 
on eqs. ( |1.1| ) might occur depending on details of the model under consideration, but 
nonetheless the equations are to a great extent universal among the various integrable 







systems. They can generally be resolved by using the Lagrange interpolation formula, 
thus leading to the formulation of a Jacobi inversion problem from which the dynamics 
of the model can be solved in terms of the zeroes of theta-functions associated with the 



curve. 



In the last few years a gradual shift of emphasis has taken place from continuous to 
discrete integrable systems. One particular class of such discrete systems represent dy- 
namical systems evolving with discrete time, namely integrable mappings, cf. e.g. || for 
a review. Such mappings can arise as periodic solutions of partial difference analogues of 
equations of KdV type, cf. |7|, ||. Since these are the discrete analogues of the finite-gap 
potentials in the continuous theory, it is natural to expect that analogous methods can 
be applied to these discrete models as the ones used to deal with periodic problems in 
the continuous case. In a recent paper Q the finite-gap integration of such mappings 
was considered and the interpretation of the maps was given as constituting special addi- 
tion formulae for hyperelliptic Abelian functions involving special winding vectors on the 
spectral curve. As a byproduct difference analogues of the Dubrovin equations ([hi]) were 
derived which form the equations of the discrete motion of the auxiliary spectrum under 
the KdV mappings. The term discrete Dubrovin equations stems from the paper (l| where 
related equations were discussed but no explicit formulae were given. 

It is well-known that in problems associated with hyperelliptic curves, such as the 
periodic problems for equations of KdV type, the poles of the Baker-Akhiezer function 



play the role of separation variables, cf. e.g. |K], 11]. Recently, this observation was 



cast into a wider framework by Sklyanin in [12] using the Lax pair approach to obtain 
separation of variables for a wide class of integrable models. In the last section of this 
note we will comment on the fine interplay between the discrete-time integrable systems 
which are the mappings of KdV type and the separation mechanism. We believe that the 
problem of finding the structure behind the discrete Dubrovin equations might form one 
of the keys to understand the separation mechanism in the new approach. 
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2 Discrete Dubrovin Equations 



In this note we investigate a class of discrete integrable systems that admit the following 
Lax description. Let us introduce the elementary matrices depending on a discrete variable 
n labelling the sites along a chain of length N, given by 

VnW = [ ^ 1 ) , (2.1) 

V K o J 

in which X n = A + a n , A being the spectral parameter and the a n being some site- 
dependent shifts. All variables can be taken to be real or complex valued as we wish. The 
most important object for us is the monodromy matrix 

JL ( A(X) B{\) \ 

T(X) =11 V n (X) = , (2.2) 

n=l V C(X) D(X) J 

which we assume to evolve under a discrete-time map v n i— ► v n according to 

f (A) = M(X)T(X)M{X)- 1 , M(A) = [ 1 | . (2.3) 

v A °y 

We will not discuss here the details of the explicit map on the level of the local variables 
v n , (see the comments in section 5), but only work on the level of the monodromy matrix. 
What is essential is that the discrete-time evolution implies the invariance under the map 
of the spectral curve 

r : $L(r], A) = det(T(A) — 77) = , n = l,...,P, (2.4) 

which defines an hyperelliptic curve of genus g = P — 1, the branch points of which being 
defined by the formula 

2<?+l 2g+l 

t = R(X) = J] TjX = r 2g+1 [] (A - e 3 ) , (2.5) 
3=0 j=i 

with 77 = 77 — |(A(A) + -D(A)). R(X) is the discriminant of the hyperelliptic curve, and its 
branch points ej, % = 1, . . . , 2g + 1 as well as their symmetric functions of order i are 
invariant under the discrete-time map Q2.3] ). In principle the branch points, which corre- 
spond to the eigenvalues of the associated N x N tridiagonal matrix, can be complex. To 



establish under what conditions they are real-valued and simple requires detailed analysis 
which is beyond the scope of this note. 

As is well-known from the continuous theory, []lO| , the roots fJLj, j = 1, ... ,g of 
the polynomial B(A), which incidentally correspond to the poles of the Baker- Akhiezer 
function, define the so-called auxiliary spectrum, and they contain the relevant information 
on the dynamics of the system under consideration. Let us, therefore, derive the discrete 
equations for the auxiliary spectrum which then form the natural analogue of Dubrovin 
equations This will, in fact, lead to a universal system of coupled first-order differ- 

ence equations associated with the hyperelliptic curve and to derive them we only need 
some global properties of the monodromy matrix. 

For genus g the monodromy matrixes T(A) takes on either one of two distinct possible 
forms depending on whether the length N of the chain in ( |2.2[) is even or odd. Writing 
eq. (|2.2|) as (in which we have taken on = in accordance with (|2.3[)) : 



T{ x , , A*+% +1 + \9 Ag + - + Ao \°B g + A^iVi + • • • + Bo ( } ( . ( 

A (\ 9 C g + \ g ~ 1 C g -i H + C ) A (\SD g + \ 9 - l D g ^ + • • • + D ) 

the two cases associated with genus g are either of the folowing two: 

n . , r / V = , B p = 1 1 / A g+1 = 1 , B g + 

Udd base : < > , bven Case : < 

{ C g = l , D g = J { C g ^0 , D g = l 

(2.7) 

The even/odd distinction is related to the variance in the choice of periodic initial data in 
the two-dimensional lattice described by lattice equations of KdV type, cf. Q, where we 
can impose periodicity on chains ( "staircases" ) with period 2P — 1 respectively period 2P 
both corresponding to a curve of genus g = P — 1. We note in passing that as a function 
of the invariants Ij, j = 0, . . . ,g of the map which are the coefficients of the trace of the 
monodromy matrix 

g 

tr T(A) = Io + Yl ^ > h = A i + Dj-l > ( 2 - 8 ) 

the top coefficient I g = A g + -D 9 -i being a Casimir with respect to the natural Poisson 
algebra associated with the dynamical map, the discriminant of the curve takes on the 



form 



R(\) = \\ lo + J2 I J xj ) -II(- A -) 

\ 3=1 J n=l 



B g C g \^ +1 + 



(2.9) 



thus depending on the shift variables a n entering via X n = A + a n . 

From the discrete-time map ( |2.3[) we have the following discrete relations for its entries 

C{\) = A-B(A) 

D{\) = A{\)-wB{\) . (2.10) 



A(X) = wB(X) + D(X) 
XB(X) = C(A) + w(A(X) - D(X)) - w 2 B{X) 



Expanding eq. ( 2.10 ) in powers of A we are lead to the following set of equations: 



A) = wB = I , Aj = wBj + Dj-i , j = 1, . . . , g 



(2.11) 



The entry B(X) has the following factorisation: 



B(X) = B 9 H(X- f i J ) , 



(2.12) 



leading to the expressions for Bo/B g , . . . , B g _i/B g as elementary symmetric functions of 
the zeroes fix, . . . ,fj, g . Similarly, the coefficients of C(A) are symmetric functions of its 
zeroes [Ax, ... ,fj, g , where the undertilde denotes the backward time-shift, and from ( 2.1C| ) 
we establish that the top coefficients B g and C g can be taken to be equal and constant. 
From the fact that 



\{A - D){X) = k\J R(X) - B(X)C(X) 



(2.13) 



where the k denotes the sign k = ± corresponding to the choice of sheet of the Riemann 
surface, subject to the condition k = —k (as follows from ( j2.13| ) and the relations ( |2.10| )), 
and taking A = [ii in ( |2.13| ) we obtain the system 

2 h - D \ ( m ■ ■ ■ 14 ^ 



/ H X -D Q \ ( 



Ml 



A ) 



(2.14) 



On the one hand, (2.14) provides us with the values of Ij — Dj-\ — Dj-\ = wBj, (j = 
1, . . . ,g), whereas on the other hand from the first of eq. ( [2.11| ) we note that w = Iq/Bq, 



hence we need the expressions for the factors Bj/Bq which can be expressed in terms 
of symmetric functions of the fij as a consequence of (§T|). Alternatively, we can take 
A = fii in (|2,13| ) and thus obtain set of first-order difference equations for the /ij, namely 
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3VT 1 - [Ky/R(n) + -I e\ +M- 1 - rZy/R(p,) - -I e ) = . (2.15) 



Mi 



In eq. (|2.15|) /x = . . . , y,gf denotes the vector with entries fij and e = (1,1,... ,1 
wheras M denotes the VanderMonde matrix 

/ 

jvc = ; : 

\ ■■■ A J 

We note that the terms with Jo can be expressed in terms of the symmetric functions of 
the /x's and ju's, namely by 



jvr 1 • 



S (-(J,! 1 ,..., -fig 1 ) , 



where S(xi, . . . ,x n ) is the vector of symmetric products of its arguments, i.e. 



Sk{ x li ■ ■ • j x g) 



E 



ll<l2<-"<lk 



Furthermore, in the spirit of ||] we can obtain from ( 2.10 ) and ( 2.13 ) by using an asymp- 
totic expansion as A — > oo a reconstruction formulae for the variable w, namely 

/ r-\ 9 29+1 

w = k(VA-Va) , A = ^2(fj.j + tij)- e i ( 2 - 16 ) 

3=1 ~ 3=1 

in terms of symmetric functions of the auxiliary spectrum. The ambiguity in the choice of 
sign reflects the reversibility of the map in the forward/backward discrete-time direction, 
and can be fixed at our choice subject to the condition that k = —k. 

The discrete Dubrovin eqs. ( [2.15 ) is a set of g first-order difference equations for 
the fij and similarly as in the continuous case (|1.1[) depending on the invariants Ij only. 
Discrete Dubrovin equations were mentioned first in 0], but not given in explicit form. 
In Q they were interpreted as special addition formulae for Abelian functions. From the 
rather general and elementary derivation given above we conclude that these equations 



ar quite universal and linked to the Jacobi inversion problem on the hyperelliptic curve 
in a fundamental way. It would be of interest, therefore, to find a mechanism similar to 
the Lagrange interpolation trick that works in the continuous case to connect the discrete 
Dubrovin equations to linear motion on the Jacobian of the curve. So far, we have not 
been able to find such a "direct" integration, although some of the original works by Abel, 
e.g. [13 1, suggest that such a mechanism might exist. We have for the time being to resort 



to an indirect approach to linearise eqs ( 2.15| ) which is by exploiting what is effectively 



the interpolating flows of the corresponding discrete-time maps. This will be discussed in 
the next section. 

We mention at this point also an alternative form of the discrete Dubrovin equations, 
namely as a second-order implicit system of equations of the form 

B(pi)-B(jii) = 2K—y/R(jii) , i = l,...,g. (2.17) 

Mi 

This form of the discrete Dubrovin equations, which is closer to the forms proposed in Q, 
was the one discussed in 01. 



3 Explicit Examples: g = 1 and g = 2 



Let us investigate the form of the discrete Dubrovin equations ( 2.15| ) in the special cases 



of genus g = 1 and g = 2. In the case of g = 1 the discrete Dubrovin equations (|2.15| ) 
reduce to one single equation, namely 



~Io (i--)= -^/W) + ^^/W) , (3-1) 



where for simplicity we have omitted the sign designation k. It is not hard to see that 
if we properly normalise the corresponding elliptic curve, i.e. ( |2.4| ) for g = 1, eq. ( [3.1| ) 
is actually nothing else than the addition formula for the Weierstrass elliptic p-function, 
namely by identifying 

V = p{8) - p(t) , Jl = p(5) - p(t + 5) , (3.2a) 

together with 

rj = = p'(t) , -~I = p'{8) , (3.2b) 
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the parameter 5 playing the role of the discrete time-step. Thus, in this case the dis- 
crete Dubrovin equation forms a natural difference equation for the Weierstrass elliptic 
p-function. Generic difference equations for elliptic functions were proposed by Potts K3] 
a decade ago. 

In the case of genus g = 2 the discrete Dubrovin equations are given by the following 
set of two coupled first-order difference equations: 



2 \MlM2 MlM2/ Ml ~ M2 

(3.3a) 



2 V /i!/x 2 M1M2 / Ml - A*2 Mi - M2 

(3.3b) 

This system can be resolved along the following lines. Without giving any details of the 
analysis we mention that the /ij arise as the poles, and the fa as the zeroes of a "transition 
factor" 92/93 where (p is the relevant Baker- Akhiezer function. In this case the independent 
variable t lives on the Jacobian Jac(T) of the curve and we have the Jacobi inversion 
problem in the form: 

/■mi a + 11b 

t = du + du , du = d/s , (3.4) 

Jei Je 2 V-R(m) 

dtt being a properly chosen vector (i.e. with appropriately chosen vectors a and b) of 
holomorphic differentials on the curve, and e\,e2 being two of the branch points. Using 
Abel's theorem the map fa — > fa is resolved via the shifts on the Jacobian of the form 

6 = 22 / du = - du e Jac(T) (3.5) 

i=l 2 f' 00 

in terms of the special winding vector (on the right hand side) on the Jacobian. 

In H we have used Kleinian functions to parametrise the solution in terms of the genus 
g = 2 analogue of the Weierstrass p-function, namely in terms of the Kleinian functions 
Pij, hj = 1)2. The definitions of the generalised Weierstrass functions go back to Klein, 



cf. |]i5| ], and have also been discussed at great length in the monographs by Baker, [16, 17 1 



The general theory has been revived in the papers [18]. (For an outline of the relevant 
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definitions, cf. the Appendix). We restrict ourselves here to list the relevant formulae for 
the case g = 2 to parametrise the solutions /x = (/ii, /i2) of the discrete Dubrovin equations 



(3.3), are the following, cf. p7| : 



P22C") = A*i + M2, Pi2(«) = -fiifJ>2, Pu{u) 



F (ni,^) - 2mrj2 
4(Mi " ^) 2 



where 



(3.6a) 



(3.6b) 



fc=0 



in terms of the coefficients of the g = 2 curve, cf. eq. fl2.4p .We note that only two of 
the functions Pij(u) are independent. The three functions P22( u ), P2i( u ) and pu(u) are 
connected by a quartic relation which is the remarkable Kummer surface % in C 3 , the 
equation of which is given by 



% : det K = det 



-ro 



2pu 



"2pl2 



\ 



in 







(3.7) 



-r 2 - 4pn ir 3 + 2pi2 2p 2 2 
2pn ^r- 3 + 2pi2 -r 4 - 4p 22 2 
\ -2p 12 2p 22 2 

The corresponding evaluations on the curve can be expressed in terms of the Kleinian 



three-indexed symbols, |L7], 16], by the relations 



Vl = A*lp222(lt) + p22l(«), V2 = A*2p222(w) + P22l(w) • 



where these odd Kleinian functions are related to the even functions via 



(3.8) 



2 

P222 



4pll + r 3 p 22 + 4p^ 2 + 4pi2p22 + r 4 p22 + T 2 , 



P221 = r - 4pnpi2 + r 4 pi 2 + 4p 22 p?2 



(3.9a) 



What is needed still is to identify the arguments of the Kleinian functions with the variable 



t on the Jacobian. This is rather subtle and will be discussed separately, cf. [19|. 



4 Connection with Separation of Variables 



As has been mentioned earlier, the auxiliary spectrum {/ij} are on the one hand the poles of 
the properly normalised Baker- Akhiezer function, on the other hand they can be considered 



to be separation variables for the integrable Hamiltonian systems, cf. [10, 11]. In recent 



papers, [12, pC|| , this connection has been put into a bigger framework to understand the 
general problem of seperation of variables for integrable systems going beyond the classical 
coordinate separation. The general statement is the following, p^| : 

For a finite- dimensional integrable Hamiltonian system ({%}, {pi}, H) admitting a Lax 
representation, i.e. for which a matrix L(X) exists such that the spectral curve 

T : det (L(A) - 77) = 

yields a complete set of integrals of the motion in involution, the poles {fa} of the Baker- 
Akhiezer function, i.e. the eigenfunction of the Lax matrix when properly normalised, 
together with the corresponding evaluations on the curve, i.e. rji = rj(fa) form a set of 
separation variables for the Hamiltonian system, i. e. there is a canonical transformation 

(Qi,Pi) ->■ (fa,Vi) 
and the separation equations are given by 

det (L(fa) - rji) = 
expressing the fact that the (fa, rji) are lying on the spectral curve. 

This statement has been exemplified by a large number of examples, cf. fl2]| , and notably 
recently by generic systems of Calogero-Moser and Ruijsenaars type even in the full elliptic 



case, [21]. We point out that the question of the normalisation of the Baker-Akhiezer 



function is a crucial point and no general prescription exists to date. 

The separation mechanism is rather independent of the type of dynamics we impose 
on the Hamiltonian system. Thus, it applies equally well to situations where we have an 
integrable discrete-time dynamics described by a Lax pair, in which case the corresponding 
discrete-time map is an (iterative) canonical transformation which usually can be cast into 
a discrete Lagrangian form, cf. e.g. H. This applies to the mappings of the type studied 



in the present paper, cf. p], and in [22] for mappings of this type a classical r-matrix 
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structure was found from which we can extract the following Poisson brackets: 

{W(A , )} _ B W AW-i WB( V) + l B(W) (41o) 

{B(A),B(V} = (4.1b) 

{W(A)} _ . M^g(W (41c) 

Eqs. ( |4.1| ) can be used to establish the canonicity of the separation variables, i.e. 

{Vi > Vj} = {lH , Mi) = , {^i , Mil = hjrii , (4.2) 

using 7](fii) = A(fii) . The symplecticity of the discrete-time map is guaranteed from the 
Poisson brackets involving the variable w: 

{A(X),w} = j{wA(X)-w 2 B(X)-C(X)) , {C(A) , w} = j (C(A) - wD(X)) 
{B(X),w} = j(A(X)-D(X)-wB(X)) , {D(X) , w} = j (C(A) - wD(X)) .(4.3) 

Thus, the general separation mechanism is consistent with the discrete-time evolution. 

Let us finish with some general remarks. It seems that the discrete-time systems and 
integrable mappings play a more profound role in the general problem of separation, cf. 
[ J23| |. In principle, the problem of finding the separation of variables in terms of discrete 
Dubrovin equations of the type ( 2.15| ), which for a iV-dimensional Hamiltonian system 



generically take the form 

G i ({Jl i },{f J , i };I 1 ,...,I N ) = , i = l,...,N (4.4) 

can be formulated as a problem of finding commuting canonical transformations. Thus, 
we might conjecture that in the relevant cases the following diagram actually represents 
a commuting diagram of integrable canonical maps as follows: 



/ s S{q l ,q i ) 



W{(JLi,(li) 
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In this diagram the S is the action functional describing the canonical transformation 
which is the discrete mapping (as can be derived from a discrete Lagrangian), and F 
denotes the generating function of the canonical transformation from the original variables 
to the separating variables. The canonical transformation, with generating function W, 
realising the dynamical mapping in terms of the separation variables, is obviously, by 
construction, an integrable map itself. Unfortunately, however, it does not seem easy in 
general to obtain explicit expressions for the generating function W, since this requires 
the elimination of the invariants entering as coefficients of the spectral curve. Noting that 
the general issue of superposition of canonical transformations and the connected problem 
of finding commuting canonical maps has to our knowledge not been addressed in full 
generality in classical mechanics, it seems that the advent of integrable discrete dynamical 
systems will make the study of these questions imminent. 



5 Conclusions 

In Q we constructed a family of exactly integrable finite-dimensional mappings from a 
lattice KdV equation by considering 'local' initial value problems on so-called 'staircases' in 
the lattice. These mappings are symplectic as consequence of the Lagrange structure of the 
original lattice KdV equation, and their complete integrability in the sense of Liouville, 
||, was established in |2^]. The resulting mappings can be written in the form of 
multidimensional rational mappings 1^ — ► l w : {v n } i— ► {v n }, n = 1, ... ,N, where the 
variables v n coincide with the ones entering in the monodromy matrix (2.2). In fact, their 



Lax description is based on the matrices V n used in section 2. In the recent paper, ||, V. 
Enolskii and the author investigated the finite-gap integration of the resulting mappings 
and their parametrisation in terms of the Kleinian functions. It was pointed out there 
that these integrable mappings have, in fact, the interpretation of being addition formulae 
for hyperelliptic Abelian functions for special winding vectors on the Riemann surface. In 
that paper we concentrated on the local description of the mappings, taking into account 
the dependence of the variables Hi on the variable n labelling the sites along the periodic 
chain as encoded in the Lax matrices V n . Interestingly, the shift variables a n play an 
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essential role in this description since they enter as distinct singularities on the curve 
defining special winding vectors corresponding to the shifts nnn + l. In the present 
paper we have not at all addressed that issue, nor the actual reconstruction formulae for 
the potential v n . We should also mention the work by Bobenko and Pinkall, [^4|], on the 
geometric aspects of lattice systems associated with equations of KdV type, where from a 
different perspective finite-gap formulae have also been presented. 
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Appendix: Kleinian Functions 



In this appendix we collect a number of relevant formulae on Kleinian functions, [15, 16, 
|j~7f . We refer the reader to the excellent review papers fTsI , (from which this material is 
taken), in which this theory has been cast in a modern context and where it was shown 
that these functions arise naturally within the KdV theory. 

Let the hyperelliptic curve be given by F : y 2 = R(x), as in ([O]), of genus g having 
2g + 2 branch points e\, . . . , e2 g +i, oo, cf. also (|2.5| ), and as is well-known we can equip it 
with a canonical homology basis (ai, . . . , a g ; b\, . . . ,b g ), and define on T a canonical set 
of holomorphic differentials, 

x^^dx 

du = (du\, . . . , du g ), duk = , (A.l) 

and differentials of the second kind with a pole at infinity 

dn = (dni, . . . , dn g ), , dUj = 22 + 1 - j>*+i+i- 4— • ( A - 2 ) 

k=j y 

Introducing g x 2g period matrices (2u,2u') and (2/7,2//) of their respective integrals 
over the a- and b cycles, we note that detu> ^ and the matrix r = uj~ 1 u' is symmetric 
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and its imaginary part is positive definite. Introducing also the matrix x = rj(2iv) , we 
can now define the fundamental hyperelliptic cr-function by the formula 

1 



a(u) 



7T9 



=exp{w'x-u}6l[e]((2^)- i u|T 



(A.3) 



where [e] 



, is the characteristic of the vector of Riemann constants, and 0[e](i?|r) 



is the Riemann theta function with characteristic, 

9[s](v\t) = ^2 exp{in(m + e') t T(m + e') + 2in(m + e'Yiv + e)} . (A.4) 

The Kleinian p-functions are defined by the second- and third order logarithmic deriva- 
tives of the Kleinian cr-function, namely 



Pij( u ) 



d 2 In a(u) 



duiduj 



, Pijk(u) 



d 3 In a(u) 



i,j,k = 1, 



,9- 



(A.5) 



duidujduk 

The Abel map 21 : (r) 9 — > Jac(T) of the symmetrised product T X • • ■ x T to the Jacobi 



variety Jac(T) = C 9 /2u © 2u/ of the curve V is defined by the vector relation 

y / du = u , 

i=l J(P,et) 



(A.6) 



and where the ej, i = 1, . . . ,g are a properly chosen selection of branch points among the 
2<7 + 2 branch points (except oo) on the curve. 

The Jacobi inversion problem of inverting the Abel map, in the formulation by Weier- 
strass, can be solved in terms of the Kleinian functions as follows: The Abel preimage of 
the point u £ JaciT) is given by the set {(yi, x±), . . . , (y g ,x g )} G (r) 9 , where {x±, . . . , x g } 
are the zeros of the polynomial 



y(x; U) = X 9 - X 9 1 pg,g(u) - X 9 ' Z p g: g-l(u) 



(A.7) 



and {yi, . . . , y g } are given by 



Vk 



du„ 



(A. 



More details can be found in [18|. 
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